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Significant decrease of spontaneous magnetization in frustrated one-dimensional ferro- and fer- 
rimagnets due to non-magnetic impurities is predicted. Using the density- matrix renormalization 
group method and the exact diagonalization method, we confirm that the total spin can vanish due 
to a single impurity in finite chains. Introducing the picture of magnetic domain inversion, we nu- 
merically investigate the impurity-density dependence of magnetization. In particular, we show that 
even with an infinitesimal density of impurities the magnetization in the ground state is reduced by 
about 40% from that of the corresponding pure system. Conditions for the materials which may 
show this anomalous impurity effect are formulated. 
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Frustrations in quantum spin systems have attracted 
much attention for its potential to exhibit new phe- 
nomena which have never been observed in unfrus- 
trated systems. Various possibilities in frustrated sys- 
tems have been suggested, such as exotic excitations 
near critical points jl|, incommensurate orders in mag- 
netic fields'^, chiral orderings0, 0, 0, |^ and disordered 
ground states /7,^]. The property we will discuss in this 
paper is also one of the phenomena where frustrations 
play an essential role, and will never be observed in un- 
frustrated systems. The property is an impurity effect 
on frustrated one-dimensional ferro- and ferrimagnets. 

Usually, a small amount of impurities has little influ- 
ence on bulk magnetic quantities, since the mean dis- 
tance between impurities is so long that the correlation 
between them is very weak and usually they affect only 
local quantities. However, in some special situations, a 
small amount of impurities can cause a bulk effect. An 
example is the impurity-induced antiferromagnetic long- 
range order (AFLRO) in quasi-one-dimensional spin- 
1/2 spin-gap systems, which was thoroughly investi- 
gated by theoretical 0, E3 numerical El| and 
experimental 0, approaches. The main feature of 
this effect is roughly explained as follows: Without an 
impurity, spins form dimers locally. By introducing non- 
magnetic impurities, moments are induced around im- 
purity sites. The moments couple one to another and 
exhibit antiferromagnetic alignment. Hence, AFLRO 
and low-energy spin- wave excitation appear in the back- 
ground of high-energy triplet excitation. For this effect, 
the spin gap and the correlations between induced mo- 
ments play an essential role. In this paper, we present 
another example that a bulk quantity, magnetization, is 
substantially influenced by a small amount of impurities 
due to a different mechanism from that of the impurity- 
induced AFLRO. The impurity effect will be realized 
even with an infinitesimal density of impurities in frus- 
trated ferro- and ferrimagnetic chains that satisfy the 
conditions we will present in this paper. Thus, in such 
systems the magnetization in the ground state will be 
significantly reduced from that of the corresponding pure 
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FIG. 1: Lattice structure of models 1 and 2 (a) without an 
impurity and (b) with an impurity at site 2i. 



systems even with a usually-negligible amount of impuri- 
ties. We notice that without the knowledge developed in 
this work, reduction of magnetization from the expected 
values tends to be explained by assuming complex, higher 
order interactions such as the Dzyaloshinsky-Moriya in- 
teraction or interchain antiferromagnetic couplings. 

In order to demonstrate the anomalous effect of non- 
magnetic impurities on frustrated ferro- and ferrimag- 
nets, we consider the minimal models defined by the fol- 
lowing Hamiltonian: 



H = Ji ^ 'S'2i-i • 'S'2i+i + J2 ^ Si • Si-^i 



(1) 



where Si denotes the spin operator at site i. The lattice 
structure is shown in Fig. Q (a). In model 1, the spin 
lengths of all spins are one half. The coupling constant 
Ji is antiferromagnetic (Ji>0), and J2 is ferromagnetic 
(J2<0). In model 2, the spin lengths at even sites are 
one half, and those at odd sites are one; both coupling 
constants are antiferromagnetic (Ji and J2>0). Model 
1 is nothing but the one proposed by Hamada and his 
coworkers in Ref. 0. Model 2 can be reduced to well- 
known models by neglecting Ji or J2: At J2=0 this model 
is equivalent to the S=l Heisenberg chain and free spins, 
and at Ji=0 it is nothing but the spin- alternating Heisen- 
berg chain. The ground states of models 1 and 2 become 
ferro- and ferrimagnetic, respectively, when IJ2I is suf- 
ficiently larger than Ji. Hereafter, the number of unit 
cells and the number of sites are denoted by L and A^g, 
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FIG. 2: (a) Schematic picture of magnetic domains (a-1) with- 
out an impurity and (a-2) with an impurity, (b) Ferrimagnetic 
chain on a bipartite lattice, where the spin lengths of all spins 
are one half, and all coupling constants are antiferromagnetic. 



respectively, and open boundary conditions are applied. 
In models 1 and 2, L=Ns/2. 

When a non-magnetic impurity is doped at an odd site, 
the system is decomposed into two pure systems. Then, 
the situation is rather trivial. Thus, we first concentrate 
on the case where an impurity occupies an even site as 
shown in Fig. ^ (b). Since the interactions from the 
impurity site are removed, the remaining interaction be- 
tween the spins adjacent to the impurity site is Ji, which 
is antiferromagnetic. Hence, we expect that the total spin 
5* in the ground state, which corresponds to the sponta- 
neous magnetization, becomes \Sl—Sr\, where Sl and 
Sr are those of the blocks to the left and right of the 
impurity site, respectively. The expected picture of this 
inversion of magnetic domains is schematically shown in 
Fig. 13 (a). As a special case where Sl=Sr^ the total 
spin will vanish due to a single impurity. 

In order to confirm this substantial decrease in spon- 
taneous magnetization, we calculated S in the ground 
states of models 1 and 2 with a single impurity put at 
all possible even sites in up to 40-site chains with even 
L by the density-matrix renormalization group (DMRG) 
method and the exact diagonalization method. The 
coupling constants are set to be Ji=0.1 and | J2| = l. We 
calculated the total spin S in the ground state by using 
the formula ^(^+1) = ((E^ (=(E^,, ^^■^J)). where 

i and j run over all sites, and ( ) denotes the expecta- 
tion value in the ground state. The numerical results on 
the total spin S satisfied the relation S=\Sl—Sr\ in all 
the cases we have investigated. There is no mathemati- 
cal proof on this relation for quantum spin systems with 
frustrations, hence it is nontrivial. The physical picture 
of the domain inversion can be intuitively understood 
by considering the corresponding Ising models, for which 
this relation holds with total spin S replaced with total 
z-component of spins . 

It should be noted, on the other hand, that in unfrus- 
trated ferrimagnetic chains on bipartite lattices such as 
shown in Fig. [21(b), this anomalous impurity effect does 
not occur, since the sign of the effective coupling between 
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FIG. 3: Impurity-density dependence of magnetization. The 
magnetizations without an impurity are denoted by Mq. Solid 
diamonds and circles denote the results on models 1 and 2, 
respectively. Open symbols denote those of the corresponding 
unfrustrated systems by setting Ji to be ferromagnetic. 



domains does not change due to impurities: In both def- 
initions of domains denoted by solid and dotted lines in 
Fig. 121(b), the sign of the effective coupling remains the 
same before or after impurity doping. Actually, in these 
systems, ferrimagnetic ground states are ensured by the 
Marshall-Lieb-Mattis theorem [T3.[T^ with or without an 
impurity. 

Based on the above numerical results for models 1 and 
2 doped with an impurity, it is natural to expect that the 
total spin S is expressed in terms of those in domains {Sk) 
as 



(2) 



when impurities are doped at even sites. Taking this 
relation into account, we have calculated magnetization 
M in an infinitesimal magnetic field with impurities ran- 
domly distributed on a chain, where impurities can sit 
not only on even sites but also on odd sites. Here, the 
magnetization M in an infinitesimal magnetic field is ex- 
pressed in terms of the total spin Si in the l-th isolated 
cluster as M=^j^Si. To be concrete, we have calcu- 
lated the average of M over 10,000 chains. Each chain 
has 100/x sites and 100 randomly distributed impurities, 
where x is the impurity density. The numerical result 
on the impurity- density dependence of magnetization is 
shown in Fig. [31 Magnetizations of models 1 and 2 are 
drastically reduced due to impurities (solid diamonds and 
circles, respectively). In particular, in the limit of small 
impurity-density, the magnetizations decrease down to 
about 57.7% of those of the corresponding pure systems; 
M(x^0):^0.577xM(x=0). 

This feature is contrasted with that without frustra- 
tions: As an example, we consider the model 1 with 
all coupling constants ferromagnetic. In this model, the 
ground state is ferromagnetic with or without an impu- 
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FIG. 4: Ground-state energies as a function of measured 
from that of 5*^=0 in a 42-site cluster of model 2 with Ji—0.1 
and J2=l without an impurity (diamonds) and with an im- 
purity at the 20-th site (circles). Solid lines are guide to eyes. 
The inset shows the calculated S. Dashed and dotted Unes 
indicate typical behaviors of a ferrimagnet and a disordered 
state, respectively. The S and of the fully polarized state 
are denoted by iSmax and 5'max? respectively. 



rity. Thus, the magnetization decreases by the amount 
of the spins at impurity sites. Namely, the magnetization 
linearly decreases as a function of the impurity density, 
i.e. M/Mo=l X ^ as shown in Fig. [31 (open diamonds). 
Since the corresponding Ising models show the same im- 
purity effect for the z-component of spins, it is expected 
that the effect shown here will also be realized in XXZ 
models with Ising-like anisotropy. 

Now, let us consider excitations in doped systems with 
non-magnetic impurities. We calculated energies of a 42- 
site chain of model 2 with open boundary conditions, 
where an impurity is put at the 20-th site. The cou- 
pling constants are set to be Ji=0.1 and J2=l. The 
finite-size algorithm of the DMRG method 17] is ap- 
plied with truncation number up to m=150. We per- 
formed 10 sweeps and confirmed convergence by calcu- 
lating 5'(5'+l) = (Ei J Si-Sj). In the inset of Fig. H we 
plot the calculated S with respect to total z-component 
of spins S^{=^- Sf). The calculated values of S almost 
coincide with the typical behaviors of a ferrimagnet and a 
disordered state (dashed and dotted lines, respectively), 
indicating that the wavefunctions with various 5'^'s are 
well converged. This figure also shows that the ground 
state changes from a ferrimagnetic state to a spin-singlet 
state due to a single impurity. The ground-state energies 
within the subspaces of fixed measured from that of 
5*^=0 are shown in Fig. 21 This figure suggests that 
there is a low-energy continuous excitation from a spin- 
singlet state 20j. This feature is contrasted with that of 
the Ising model where the lowest excitation has a finite 
gap of the order of Ji or J2 independently of the clus- 
ter size. The ground states with small S'^'s in doped 
quantum systems are almost degenerate, which would be 
reflecting ferromagnetic fluctuations in domains. 



FIG. 5: Magnetization curve of model 2. The parameters are 
the same as those in Fig. 0] 



We calculated magnetic field H by using a discretized 
form of the derivative of energy E with respect to S^: H= 
dE/dS^c:,{EiS^^,)-E{S:,)}/{S^^+,-S^,}, where S^=n 
or n+0.5 with or without an impurity. (n=0, 1, • • • .) 
The result on the magnetization curve is shown in Fig. 
ISl The magnetic field required for the magnetization to 
recover up to the spontaneous magnetization of the pure 
system is about 0.1 which is the order of Ji as expected 
from the picture of domain inversion (Fig. [21(a)). 

Based on the above considerations, we list the condi- 
tions for the anomalous impurity effect: 

1. The system should be one-dimensional. Namely, 
interactions between chains should be much smaller 
than those in chains. 

2. The ground state without an impurity should have 
spontaneous magnetization. 

3. Local interactions near impurity sites should be set 
such that the effective interaction between mag- 
netic domains changes from ferromagnetic to an- 
tiferromagnetic due to impurities. 

The third condition leads to frustration. 

The models that satisfy the above conditions will ex- 
hibit the anomalous impurity effect. For example, the 
decorated triangle chains (Figs. IHl (a-1) and (a-2)) and 
the diamond-like chain (Fig. ^ (b)) will be the models 
that exhibit this effect with all spins one half and all 
coupling constants antiferromagnetic. In the decorated 
triangle chains, when J2 is sufficiently larger than Ji, 
the spins on decorating sites align parallel, resulting in 
a ferrimagnetic ground state. If an impurity is doped on 
the top site of a triangle, the remaining interaction be- 
tween the spins adjacent to the impurity is Ji, which is 
antiferromagnetic. Thus, the domain inversion and sub- 
stantial decrease in magnetization due to impurities are 
expected. Actually, we have confirmed by exact diagonal- 
ization that the total spin S in the ground state behaves 
as S=\Sl—Sr\^ when an impurity is doped on top sites 
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FIG. 6: Possible models for the anomalous impurity effect 
with spins one half and coupling constants antiferromagnetic. 
(a) Decorated triangle chains, (b) Diamond-like chain. 

of triangles in up to 24-site clusters with Ji=0.1, J2 = 1.0 
and J3=0.5|2l|. 

In the case of the model in Fig. EKb), the ground state 
becomes ferrimagnetic, when J2 is sufficiently larger than 
Ji. If an impurity is doped at site 3z, the effective inter- 
action between the spins at sites 3i—2 and 3i-\-l is mainly 
determined by the three-site Hamiltonian of sites 3i— 2, 
3i—l and 3z+l. In order for the effective interaction to be 
antiferromagnetic, Ji has to be larger than the effective 
coupling by J2S through the spin at site 3i — 1. If such a 
parameter can be chosen, the magnetic domain inversion 
due to impurities will be realized. In this paper, we do 
not intend to determine the precise boundaries for this 



effect, since in delicate systems such as that of Fig. ini(b) 
the phase boundary will depend on the system size. In- 
stead, we would like to emphasize that, as demonstrated 
in this paper, there actually exist systems that exhibit 
this impurity effect in some parameter regimes for frus- 
trated ferro- and ferrimagnets in one dimension. 

In summary, we have investigated effects of non- 
magnetic impurities on frustrated ferro- and ferrimagnets 
in one dimension by the DMRG method and the exact 
diagonalization method. Based on the numerical results, 
we pointed out that in these systems a small amount of 
impurities can drastically decrease magnetization in the 
ground state. Introducing the picture of magnetic do- 
main inversion, we have investigated impurity-density de- 
pendence of magnetization. In particular, we have shown 
that the magnetization with an infinitesimal density of 
impurities becomes as small as 57.7% of that without an 
impurity. The energy scale of this impurity effect is of 
the order of the remaining effective interaction between 
the spins adjacent to impurity sites. The low-energy exci- 
tations in doped systems are continuous from the lowest 
spin-state (except the finite-size gap). We also listed the 
conditions for this impurity effect. 

In the materials which are effectively described by 
frustrated spin models, other interactions such as the 
Dzyaloshinsky-Moriya interaction or biquadratic interac- 
tions are sometimes not negligible. Although their influ- 
ence on the impurity effect requires further study, the 
prediction in this paper deserves careful experimental in- 
vestigations. 

We thank M. Hase for valuable discussions and com- 
ments on related materials. 



[1] T. Senthil, A. Vishwanath, L. Balents, S. Sachdev and 

M.P.A. Fisher, Science 303, 1490 (2004). 
[2] N. Maeshima, K. Okunishi, K. Okamoto and T. Sakai, 

Phys. Rev. Lett. 93, 127203 (2004). 
[3] D. Grohol, K. Matan, J.-H. Cho, S.-H. Lee, J. W. Lynn, 

D. G. Nocera and Y.S. Lee, Nature Mat. 4, 323 (2005). 
[4] S. Fujiki and D.D. Betts, Can. J. Phys. 65, 76 (1987); 

Prog. Theor. Phys. Suppl. No. 87, 268 (1986). 
[5] H. Nishimori and H. Nakanishi, J. Phys. Soc. Jpn. 57, 

626 (1988). 

[6] P.W. Leung and K.J. Runge, Phys. Rev. B 47, 5861 
(1993). 

[7] P.W. Anderson, Mater. Res. Bull. 8, 153 (1973). 

[8] P. Fazekas and P. Anderson, Philos. Mag. 30, 432 (1974). 

[9] H. Fukuyama, T. Tanimoto and M. Saito, J. Phys. Soc. 

Jpn. 65, 1182 (1996); H. Fukuyama, N. Nagaosa, M. 

Saito and T. Tanimoto, ibid. 65, 2377 (1996). 
[10] M. Sigrist and A. Furusaki, J. Phys. Soc. Jpn. 65 2385 

(1996). 

[11] N. Nagaosa, A. Furusaki, M. Sigrist and H. Fukuyama, 

J. Phys. Soc. Jpn. 65, 3724 (1996). 
[12] Y. Motome, N. Katoh, N. Furukawa and M. Imada, J. 

Phys. Soc. Jpn. 65, 1949 (1996). 



[13] Y. lino and M. Imada, J. Phys. Soc. Jpn. 65, 3728 (1996). 
[14] M. Hase, K. Uchinokura, R.J. Birgeneau, K. Hirota and 

G. Shirane, J. Phys. Soc. Jpn. 65, 1392 (1996); M.C. 

Martin, M. Hase, K. Hirota, G. Shirane, Y. Sasago, N. 

Koide and K. Uchinokura, Phys. Rev. B 56, 3173 (1997). 
[15] M. Azuma, Y. Fujishiro, M. Takano, M. Nohara and H. 

Takagi, Phys. Rev. B 55, R8658 (1997); M. Azuma, M. 

Takano and R.S. Eccleston, J. Phys. Soc. Jpn. 67, 740 

(1998). 

[16] T. Hamada, J. Kane, S. Nakagawa and Y. Natsume, J. 
Phys. Soc. Jpn. 57, 1891 (1988). 

[17] S.R. White, Phys. Rev. Lett. 69, 2863 (1992); Phys. Rev. 
B 48 10345 (1993). 

[18] W. Marshall, Proc. Roy. Soc. A 232, 48 (1955). 

[19] E.H. Lieb and D. Mattis, J. Math. Phys. 3, 749 (1962). 

[20] In finite-size systems, there should be a gap due to the 
discreteness of energy levels. Here, we mean that an 
asymptotically continuous low-energy excitation exists 
which reflects the quantum nature of the system. 

[21] For the model in Fig. (a-2), the degrees of freedom of 
the spin on the decorating site connected to the impurity 
are neglected, since it behaves as a free spin. 



